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A thermal relic, often referred to as a weakly interacting massive particle (WIMP), is a particle 
produced during the early evolution of the Universe whose present (relic) abundance depends only 
on its mass and its thermally averaged annihilation cross section (annihilation rate factor) {(Tu)ann- 
Late time WIMP annihilation has the potential to affect the cosmic microwave background (CMB) 
power spectrum. Current observational constraints on the absence of such effects provide bounds on 
the mass and the annihilation cross section of relic particles that may, but need not be dark matter 
candidates. For a WIMP that is a dark matter candidate, the CMB constraint sets an upper bound 
to the annihilation cross section, leading to a lower bound to its mass that depends on whether or 
not the WIMP is its own antiparticle. For a self-conjugate WIMP, rrimin = 50 / GeV, where / < 1 
is an electromagnetic energy efficiency factor. For a non self-conjugate WIMP, the minimum mass 
is a factor of two larger. For a WIMP that is a subdominant component of the dark matter density 
there is no bound on its mass and the upper bound to its annihilation cross section imposed by the 
CMB transforms into a lower bound to its annihilation cross section. These results are outlined and 
quantified here using the latest CMB constraints for a stable, symmetric (equal number of particles 
and antiparticles), WIMP whose annihilation is s-wave dominated, and for particles that are, or are 
not, their own antiparticle. 


I. INTRODUCTION 

Any new, massive, beyond the standard model particle, Xj populated by interactions with the standard model 
particles present during the early evolution of the Universe and in thermal equilibrium with them, will annihilate when 
it becomes non-relativistic (for temperatures T < m^), reducing its abundance (relative to, e.g., photons). Initially, 
the annihilation rate (per particle) Fann = n^{o'v)i,nn cc T^, exceeds the expansion rate, the Hubble parameter, 
H (X and = n^^eq, the number density in thermal equilibrium at temperature T. Annihilations continue 
until the annihilation rate (per particle) becomes too slow compared with the expansion rate of the Universe, and 
the abundance “freezes out”, when T = Tf, where Fann,/ = Hf- For T < Tf, Fann < H, preserving the relic 
number of particles in a comoving volume. For symmetric {n^ = n^), stable particles, the late time abundance {e.g., 
the ratio of the current mass density to the present critical mass density, = (Px/Pcrit)o) is determined by the 
annihilation cross section, {o'v)y^ = ((Tu)ann Such particles are referred to here as thermal relics or, interchangeably, 
as weakly interacting massive particles (WIMPs). In the WIMP “paradigm”, the annihilation cross section is chosen 
{{o'v)y^ = (err’) dm) so that the WIMP accounts for the observationally inferred dark matter (DM) density, = Udm- 

However, annihilations don’t cease when T = Tf. It is important to realize that even though when T < Tf, 
Fann < H, the WIMP continues to annihilate throughout the evolution of the Universe, independent of whether 
or not it is a dark matter candidate. Although long after freeze out (when T ^ Tf ^ ’^x) residual annihilations 
are rare (Fann ^ H), the energy injected from each annihilation, 2m^, may be large compared to the thermal 
energy, leading to observational constraints on the late time annihilations. Some of the earliest constraints on late 
time annihilation have come from big bang nucleosynthesis (BBN) [1-4]. The BBN constraints were supplemented 
(superseded) in a series of papers [5-8] using constraints on the deviation from black body of the spectrum of the 
cosmic microwave background radiation (CMB) frequency spectrum [9] and from the reionization and heating of the 
intergalactic gas [10]. For a discussion and references to complementary constraints from laboratory experiments, the 
effect on stellar evolution, from gamma ray observations of the Galaxy and beyond, see, e.g., [8]. At present, the 
Planck satellite observations of the CMB power spectrum [11] appears to offer the most restrictive constraints on late 
time annihilations and, in particular, on the annihilation cross section, (cru)ann [12-17]. For a dark matter candidate 
(z.e., if = Ddm)i the Planck CMB observations, by requiring that (cru)ann = (o’'i’)x = (o'^^)dm < (o’u)cmb cx m^, 
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set a lower bound to the x mass, mmin, in the interesting range of a few tens of GeV (see, e.g., [17]). The most recent 
Planck results [18] improve on these results, increasing the value of TOmin by more than a factor of two. 

In this paper these constraints from the CMB on and {crv}-^ = ((Tu)ann nre revisited, noting that there are 
different constraints depending on whether the thermal relic particle is, or is not, its own antiparticle (x = X or 
X ^ x) ) investigating how the constraints change if the WIMP is only a subdominant component of the dark matter 
< IIdm)- In particular, it is noted here that if < SIdm, then the lower bound on the WIMP mass from the CMB 
disappears, and the upper bound on (cru)ann = (o'^)x becomes a lower bound, (uv)^ > (cru)min = ((o'I')dm)^/(c’'^^)cmBj 
where (err !)dm is the value of the annihilation cross section, (cru)^, required for = Udm- An overview of the 
calculation of the thermal relic abundance and its connection with the particle mass and annihilation cross section 
is presented in §11, as a prelude to considering in §III the relation between the annihilation cross section and 

the CMB cross section (cru) cmb that emerges from the CMB constraint on late time annihilations at the epoch of 
recombination. The results presented here are summarized and discussed in § IV. 


II. THE THERMAL RELIC ANNIHILATION CROSS SECTION AND THE RELIC ABUNDANCE 

In this paper, a WIMP is a particle whose relic abundance is determined by its thermally averaged annihilation 
cross section (annihilation rate factor, (cru) ann)- This excludes from consideration here any particle - antiparticle 
asymmetry (n^ ^ n^), since in this case the relic abundance is determined by the adopted asymmetry {e.g., by 
the chemical potential) and not directly by the annihilation rate factor. Also excluded from consideration here is a 
long-lived, unstable particle whose late time energy injection depends on its mass and lifetime. This latter case has 
been analyzed in [19]. 

For a thermal relic, during the early evolution of the Universe when the temperature exceeds the WIMP mass, 
annihilations are compensated by the inverse processes (particle creation), and the WIMP abundance is preserved as 
the Universe expands and cools. In this early, extremely relativistic (ER) regime, the ratio of x particles to photons, 
n^/n^ = n^^eq/n-y = constant since, in this ER regime, n^,eq oc oc n^. When T ^ m^, particle creation is 
suppressed by an energy barrier and the thermal relic abundance decreases exponentially {n^^gqlrij oc , where 

X = m^/T). Nonetheless, depending on the strength of the coupling between the x and the standard model particles 
whose masses are < T, in this intermediate non-relativistic (NR) regime, Pann H , and the thermal relic abundance 
remains close to equilibrium, decreasing exponentially as the Universe expands and cools. This continues 

until the inverse reactions become too slow to maintain equilibrium (z.e., when the deviation from equilibrium, 
(n^ — n^^eq)/nx^eq, grows to of order unity, at T = when {n^^eq/n-y)* oc <C 1). Thereafter, for T <T^, 

annihilations dominate, quickly reducing, even further, the abundance of the x particles, until the annihilations finally 
“freeze out” at T = Ty « T^/2, when Tann,/ = (’^x(®’^)ann)/ = Hf. For a thermal relic, the WIMP abundance {e.g., 
the ratio of the current WIMP mass density to the present critical mass density, = (Px/Pcrit)o) is determined by 
the annihilation cross section, (cru);,^ = (crz;)ann- 

The time (temperature) evolution of a symmetric (n^ = n^), stable, thermal relic follows the evolution equation 
first described by Zeldovich [20], 

dN^/dt = dN^/dt = (cru)ann {nx,eq “ ^x)^ = (o■^')ann (’^x.e? “ ’ (1) 

where = n^V is the number of x particles in the comoving volume V at time t and is their number density, 

and similarly for the antiparticle y. For the analysis presented here, attention is focused on a relic particle whose 
annihilation is s-wave dominated^. For this choice, the annihilation rate factor (“cross section”), ((Tu)ann = (o’^)x! 
is independent of temperature or velocity, resulting in a cross section that is the same at freeze out and during all 
subsequent epochs in the evolution of the Universe as well as in all environments in the present Universe {e.g., the 
Sun, galaxies, clusters of galaxies, etc.). The thermal relic evolution equation, a form of the Riccati equation, cannot 
be integrated in closed form, but needs to be solved either numerically or by using an analytic approximation first 
outlined by Zeldovich [20] . The quantity that is actually predicted by the numerical and semi-analytic solutions of the 
evolution equation is the product of three terms: the mass, the annihilation cross section, and the frozen out abundance 
of the relic particles. The latter may be written in terms of the present {z = 0) number density n^o = = n^, so 

that post freeze out, for z < Zf, where Zf corresponds to T = Tf, n.^{z) = n^o{l + z)^. The solution to the evolution 


^ At late times and low temperatures {e.g., around recombination), p-wave annihilation is suppressed by a factor of Tjra^^ (or (u/tiy)^) 
that is typically <g 1, leading to much weaker CMB constraints for this case. “Sommerfeld” or resonance enhancement (or suppression) 
of the annihilation cross section is model dependent and is ignored in the discussion here. 
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FIG. 1: (Color online) The s-wave, thermal annihilation cross section (rate factor), ((tw)dm = (c’'n)ann = {(7v)x, as a function of 
the relic particle mass, required for the WIMP to account for the observed dark matter density = floMh^ = 0.12). 

The top (purple) curve is for a non self-conjugate “NS” particle (x 7^ x) and the bottom (blue) curve is for a self-conjugate 
“S” particle (x = x)- Only along the curves does the WIMP account for the DM (D^ = Ddm). Although the regions above 
the curves are allowed, in these regions the WIMP is a subdominant contributor to the DM (D^ < ^^dm). The regions below 
the curves are forbidden since, in these regions the relic WIMP mass density exceeds the DM density (D^ > Ddm). 


equation predicts the combination Airrix) = ™x^xo('^^)x [sse, e.^., [20] or [21]), where A{mx) is a slowly varying 
function of the mass (see, e.g., [22, 23]). 

The quantitative relation between the annihilation cross section and the relic mass density depends on whether 
or not the particle is its own antiparticle (x = X or x 7^ x)- When x = X) th® particle will be referred to as 
“self-conjugate” (S), while for the case where x ^ X tti® particle will be referred to as “non self-conjugate” (NS). 
To fix the normalization for the (cru)^ — f^x^^ relation^, a choice needs to be made. For S particles the mass and 
number densities are related by p^s = ''tTx^xSj while for NS particles p^NS = ™x(^xNS + ''^yNs) = 2mxnxNS- As a 
result, there are two possible choices. For example, it might be useful to choose to have the relic number densities 
equal, independent of whether or not the particle is its own antiparticle. In this case, for n^s = the solution 

of the evolution equation (evaluated at the same mass) requires that (ctu)n = (ctu)ns- For this choice the relic mass 
densities differ, p^NS = 2pxS- However, the more common choice is that, independent of whether the relic particle is 
or is not its own antiparticle, the relic mass densities are the same, p^NS = PxS- In this case, ((Tu)ns = 2 (cru)s, and 
the abundances (ie., the post freeze out number densities) of the thermal relics differ, n^s = 2nxNS- 

For the choice made here, Hns = ^Si the results for {crv)^ = (o'u)dm are shown by the two curves (for S and NS 
relic particles) in Figure 1 for fl^h^ = HDMh^ = 0.12 [11]. As may be seen in Figure 1, for a fixed value of the present 
mass density contributed by a thermal relic (x), the required cross section is not, independent of the particle mass, 
as is often assumed/claimed. The results shown here (and earlier in reference [22]) are in quantitative agreement with 
those that follow from the publicly available DarkSUSY code [24] (P. Gondolo, Private Communication). 

After freeze out (z < z/), the number density of relic particles varies with the redshift as (H-z)^, as does the number 
density of baryons (the relic of an asymmetry between baryons and antibaryons in the early Universe), so that for 
z < Zf, rixIn-B = constant. Post freeze out, the annihilation rate per relic particle is Fann = ’^x('^'^)x “ 
so that when evaluated at the same mass, Fann ,ns = Fann.s- However, the late time annihilation rates per unit volume 


® It is convenient (and conventional) to parameterize the present mass density by D^h^, where = Px/Pcrit, Pcrit is the present Universe 
critical mass density, and the present value of the Hubble parameter is Hq = 100hkms“^ Mpc“^. 






4 


differ between S and NS thermal relics. The annihilation rate per unit volume is n^^Tann = ny-A(jn^)/m^, so that, 

(^X^ann)NS _ ^^xNS _ 1 / PxNS ^ ^ 

ann )s n^s 2 V PxS / 


For the choice adopted here, PxNS = PxSj requiring ((Tu)ns = 2((Tu)s, this leads to (np^rann)s = 2(n^rann)NS- However, 
since in each annihilation, an energy E = 2m^ is injected (the same for S and NS), the energy injection rate per unit 
volume, de/dt = 2m^{n^r^nn)- But, provided that = PxSi {de/dt)j<!S = PxNsBann.NS = PxS^ann,s = {de/dt)s. 
Thus, although the late time annihilation rates per unit volume (or, per baryon) differ between S and NS particles, 
the energy injection rates per baryon are the same for both cases. 

The results for the dependence of the annihilation cross sections on the WIMP mass shown in Figure 1, spanning the 
mass range 100 MeV < < lOTeV^, are for the case where the WIMP accounts for the dark matter (fl^ = Hdmj 

where = 0.12 [11]). Only those WIMPs whose s-wave annihilation cross section and mass lie on one or the 

other of the two curves in Fig. 1 can account for the observationally inferred dark matter mass density. These curves 
should have a finite thickness, not shown in the figure, to account for the theoretical uncertainty in the (ctvJdm ~ 
WIMP mass relation (~ 5— 10% [22]) and the uncertainty (smaller) in ("^ 1.5 — 2% [18]). As may be seen from 

Fig. 1, for either kind of relic particle, S or NS, (cru)DM is very nearly constant (independent of mass) for > 10 GeV, 
with the ratio of cross section values at the same mass scaling with the factor of two explained above. However, for 
< 10 GeV, the required annihilation cross section, (crw)DM) increases by a factor of two or more, a reflection of 
the changing number of degrees of freedom contributing to the energy and entropy densities as the relic particle is 
freezing out (see, e.ff., the discussion in [22]). As noted above, the quantity predicted by the numerical and semi- 
analytic solutions of the evolution equation that lead to the results shown in Fig. 1, is the product of the annihilation 
rate factor (crv)^ and the mass density parameter, H^h^, (see, e.ff., [20] or [21]). In general, for ^ Hdm (but for 
not too different from UdmIi^), Wx(^xb2) « (o'u)DM(HDMh^)- When ^x^ differs signihcantly from 
there is a logarithmic correction to this approximation [22], 


f Mx ^ f »xh^ \ « 1 _ liog f« 1 -p ilog f 


(3) 


For values of (crt’)x and along the curves in Fig. 1 (((Tt')ann = (cru),^ = (ctvJdm), ^x ~ ^dm, while the regions 
above (below) the curves correspond to, < flnm (^x > ^dm)- The regions below the curves (for each kind of 
particle) are excluded, since in these regions > Hdm- The regions above the curves are allowed, but the particles 
corresponding to these combinations of cross section and mass will not account for the dark matter, but will only be 
subdominant contributors to the total dark matter mass density (H,^ < Hdm)- 


III. THE CMB CONSTRAINT ON ENERGY INJECTION FROM LATE TIME ANNIHILATIONS 

After the WIMP has frozen out at redshift Zf, its number density redshifts, n^{z) = (l-l-z)^n^o, so that at redshift 
z < Zf the annihilation rate per unit volume is n^{z)'^{av)^ = (1 -I- z)^n^Q{av).^. Since the relic number density may 
be written in terms of the relic mass density and the mass, and the relic mass density may be written in terms of 
the mass density parameter the annihilation rate per unit volume at redshift z < Zf is proportional to the 

combination (1 -P z)^{il.^h.^/m^)'^{av)^. Assuming that none of the annihilation energy goes into neutrinos or the 
dark sector, the effect on the background plasma of late time annihilations at redshift z due to the electromagnetic 
energy injection rate per unit volume is proportional to 

de^/dt oc (1 -k zf{n^h^ff{z){av)^/m^ , (4) 

where f{z) is an electromagnetic energy efficiency factor. This expression is general in the sense that there is no 
assumption here that the WIMP is the dark matter. 

For particles in the mass range considered here, the relic abundance is frozen out at temperatures above a few 
MeV, well before recombination, while attention here is focussed on the energy injected at the recombination epoch 
(zrec ~ 1100), where Tree ~ 0.3eV <C m^. Because the electromagnetic energy injected during this epoch has the 


^ Here, masses and mass densities are measured in energy units (c = 1). 
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FIG. 2: (Color online) The annihilation rate factors {crv)x = (o'w)dm from Fig. 1 for a WIMP that is the DM are shown 
as functions of the WIMP mass, along with the CMB constraints on late time annihilations, (o'w)cmb (shown in black) for 
/ = 1 (solid) and / = 0.2 (dashed). The energy injection from late time annihilations is constrained by (cru)cMB in that for 
ilx ~ Ddm, the total annihilation cross section satisfies (cru)^ = ((ju)dm < (crujcMB, leading to a lower bound to the WIMP 
mass where the curves cross. But, for < Ddm, there is no constraint on the WIMP mass and (cru)^ < ((tw)cmb(Dj<./Ddm)^. 
See § III A for details. 


potential to affect the CMB power spectrum [12-15, 17, 18], the observed CMB power spectrum sets a constraint on 
late time annihilations, leading to a constraint on the combination {ily^h'^)‘^f{z){{av)-^/m^) [12-15, 17, 18]. 

Since the holy grail of particle cosmology is to identify a particle that is a potential dark matter candidate, in the 
literature it is almost always assumed that the WIMP is the DM (12^ = SIdm), so that f{z){{av)^/m^) = 

f{z){av)Y)M/nT'x- This latter combination is usually identified as Pann = /(- 2 )((Tu)dm/w;^, and the CMB observations 
are used to set an upper bound to Pann- It is convenient to introduce a CMB cross section (rate factor) related to Pann 
by, ((Tr')cMB = rn^pannl f{z). For the current Planck [18] results, (cru)cMB = 4 x 10 “^®TOj^// (in units of cm^s“^, 
with in GeV). This cross section is shown as a function of by the black curves in Figure 2 for two choices of 
/. Note that while in general the efficiency factor is redshift dependent, the quantity that enters here is an “effective” 
efficiency factor, / « f{z = 600) (see [14, 17] for discussion and further references). According to [17], it is likely that 
0-2 ^ ^ 1 and these two choices are shown by the solid (/ = 1) and dashed (/ = 0.2) curves in Fig. 2. 

A naive, albeit incorrect interpretation of Fig. 2 would be that only the wedge-shaped regions above the (cru) dm 
curves and below the (ctvJcmb curves are allowed for consistency with the mass density constraint, < fluM, and 
with the CMB constraint on late time annihilations. Indeed, for the case where the WIMP is the DM (D^^ = Ddm), 
there is a lower limit to the WIMP mass, identified in Fig. 2 by the values of the masses at the crossings of the 
(crv)^ and the (cru)cMB curves. For a self-conjugate WIMP, the minimum mass is mmin ~ 50 (10) GeV for / = 1 (0.2), 
and the minimum masses are twice as large for a non self-conjugate WIMP. A stable, symmetric WIMP that accounts 
for the DM, whose annihilation is s-wave dominated, must have > rumin if it is to be consistent with the CMB. If 
future CMB observations should reduce the current upper bound on (crw)cMB [18], the lower bound on the minimum 
mass of a dark matter candidate will increase. However, as explained below in §IIIA, if the WIMP does not account 
for all of the DM and < Hdm is allowed, there is no restriction on the WIMP mass and the only restriction on 
the annihilation cross section is that it exceed some minimum value ((av)^ > (crv) min)- 
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FIG. 3: (Color online) The minimum annihilation rate factor consistent with the CMB constraint for / = 1 and the mass 
density constraint is shown as a function of the WIMP mass. The regions above the curves, corresponding to 

< Udm, are allowed, while the regions below the curves, corresponding to > Udm, are excluded. The upper (purple) 
curve is for an NS WIMP (x 7 ^ x) the lower (blue) curve is for an S WIMP (x = x)- For ^ rumin, (o-v)min = (cru)DM, 
while for < rUmin, (uu)niin = ((<^v)dm)^/(o'v)cmb > {uw)dm, corresponding to < Udm- See the text for details. 


A. Constraints On A WIMP That Is Not The Dark Matter (11^,. < Udm) 


As seen in Fig. 2, if the WIMP is the DM, the CMB constraint on Pann, requiring that (av)^ = (fTu)DM < (o'u)cmb, 
leads to a lower bound to the WIMP mass. However, allowing for a WIMP that is only a subdominant contributor 
to the DM mass density < Hdm); the CMB constraint on pann leads to a constraint on the WIMP annihilation 
cross section that not only depends on (cru)cMB, but also on the ratio 

/ \ ^/ \ 

(av)^ < (crv) CMB H h 2 7 ■ 

The CMB constraint does not limit the mass of a WIMP whose contribution to the DM is subdominant, but it 
does set a lower bound to the annihilation cross section for such a particle. Up to the logarithmic corrections in Eq. 3 
(included later), the upper bound to the annihilation cross section from the CMB, translates into a lower bound to 
the annihilation cross section, as may be seen by combining Eqs.3&5, 


(cru)x > (crvjmin = 


((crvjpM)^ 

(^v)cMB 


or 


{<^v)x > (o'u)niin 

(c^^')DM “ (o-u)dM 


(o'u)cMB 


( 6 ) 


For > TOinin, all values of (cru);^ > (cru)DM are allowed since they correspond to < Ddm, and they also satisfy 
the CMB constraint (Eq. 5)®. Another consequence of allowing < 0 dm is that the lower bound on the WIMP 
mass disappears. For all < rrimin, (<7u)„iin = ((o'u)dm)^/(o’u)cmb, and all values of (uv)^ > (CTu)„iin are allowed 
since they, too, are consistent with the CMB constraint (Eq. 6). 


® For > mmin, the CMB constraint may be rewritten as (f7n}^(fIj<.h^/r2DMh^)^ ~ (cr-i;}DM(IIxF^/^DMh^) ^ {crnlDM < {'rv}cMB (see 
Fig. 2 in support of this last inequality for > mmin)- 
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FIG. 4: (Color online) The two panels show the results for the minimum annihilation rate factor for a self-conjugate (S) thermal 
relic, (crujmin, satisfyin g the CMB constraint and < floivi, as in Fig. 3. In both panels the dotted curves account for the 
logarithmic corrections to ((Tu)min discussed in the text. The regions above the solid (dotted) curves are allowed, while the 
regions below them are excluded. In the left hand panel the solid (dotted) curve is for all masses and for / = 1. The dashed 
curve shows the extension of (nujDM to lower masses, < mmin, illustrating that (av)min > (o'v}dm in this mass range. The 
right hand panel compares ((Tw)min for f — 1 (solid) and / = 0.2 (dashed). 


These results are shown in Figure 3 for S and NS thermal relics, where (crvjmin (see Eq. 6) is shown as a function of 
(for / = 1). For an NS particle, for all masses > TOmin = 100GeV, (av}min = (o’u)dm « 4x 10“^® cm® s“^. For 
masses below this value, (crvjmin increases with decreasing WIMP mass approximately as as seen in Figs. 1 & 2, 
modulo the variation of (ctvJdm with mass shown in Fig. 1 and the logarithmic corrections from Eq. 3 described below. 
For an S particle, the same behavior for (crw)niin as a function of seen for an NS particle is shifted in mass and 
normalization. For an S particle, mmin = 50 GeV, and for > 50 GeV, ((J?;)min = (o'u)dm « 2 x 10“^®cm®s“®. 
Note that when > 100 GeV, (cru)niin,Ns/(o’t!)niin,s = 2. In the mass range from 50 to 100 GeV, (cr'c)min.NS 
increases approximately as m~^, while ((Tu)inin,s is approximately constant. As a result, when = 50 GeV, 
(o'u)niin,Ns/(o’u)niin,s ~ 4. Thereafter, for < 50 GeV, the ratio ((Tr!)„iin,Ns/(o'i')min,s = 4, remains constant for all 
lower masses. As confirmed by Fig. 3, when allowing for < fiuM, there is no bound on the relic particle mass, 
although there is a lower bound to the annihilation cross section. 

The results outlined here and shown in Fig. 3 are expanded upon in the two panels of Figure 4 where, for an 
S particle, the corresponding (avjmin curves are shown as functions of the WIMP mass, including the logarithmic 
corrections from Eq. 3 that were ignored in Fig. 3. In the left hand panel, (cru)niin and (ctvJbm are compared for / = 1. 
In the right hand panel, the results for (cru) min are shown for / = 0.2 as well as for / = 1. 

Including the previously neglected logarithmic corrections results in 


(crv}r 


((o'v)dm) 

(o'u)cMB 


2 r 


1,1, / (o-z;)dm 

1 + o log 7 — j - 

8 \(a-u)cMB 


< 


((q'ti)DM)^ 

(o'r')cMB 


or 


(cru)r 


(av) 


DM 


(crv) 


DM 


(av) 


CMB 


1,1, (^v)dm 
1 -I- - log ' 


(crv) 


CMB 


< 


(av) 


DM 


(av) 


CMB 


(7) 

( 8 ) 


For / = 1 the ratio of cross sections, (cr'(;)min/( o'i^)dm, is shown in the left hand panel of Figure 5. The regions above 
the curves ((cr?;)^ > (cru)min) are allowed, consistent with < Udm and with the CMB constraint on late time 
annihilations. 
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FIG. 5: (Color online) The left hand panel shows (for / = 1) the ratio of the minimnm annihilation rate factor to the DM 
annihilation rate factor, (crv)min/{o'n)DM, as a function of the WIMP mass. The regions above the curves, corresponding to 
> (crv}^in > (o'n)DM, are allowed, while the regions below the curves are excluded. The right hand panel shows the 
variation of the ratio of mass densities, flmax/DDM, with WIMP mass (also for / = 1). For the ratio of mass densities, the 
allowed regions are restricted to lie below the curves (fl^ < Umax < flDiu). In both panels the dotted curves account for the 
logarithmic corrections to ((Tw)inin discussed in the text. 


When (crv)^ > (cryjaiin, there is an upper limit to the relic mass density, < Umax < ^^dm, where 

^max ^ (o't;)cMB ^ 1 ^ (o't;)DM 

(o'u)dM L 8 V('^'*^)cMB 

This ratio is shown for an S particle, for / = I, as a function of the relic particle mass in the right hand panel of 
Figure 5. Here, the regions below the curves (fl^ < flma.x) are allowed. 

IV. SUMMARY AND DISCUSSION 

Although a thermal relic particle freezes out during the early evolution of the Universe, it doesn’t cease annihilating; 
“annihilations are forever” [3]. At freeze out, when T = T/, the annihilation rate (per particle) is equal to the expansion 
rate, Fann,/ = Hf. Throughout the subsequent evolution of the Universe, for T <Tf, Fann < H, and annihilations 
do continue, albeit at an ever decreasing rate (Fann ^ H). However, since the energy released in each annihilation, 
2m^, can be very large compared to the thermal energy (the temperature, T < Tf ^x) universal plasma, 

the energy released by the electromagnetic component of these rare, late time WIMP annihilations may leave an 
imprint on the CMB frequency or power spectrum [8, 12-15, 17]. Indeed, the current, best constraints are inferred 
from observations of the CMB power spectrum [18]. The CMB constraint, ((Tu)cmb < 4 x f cm.^ s~^, when 

compared to the WIMP annihilation cross section required if the WIMP is a DM candidate (U^^ = Udm), results in 
setting a lower bound to the mass of the DM WIMP (see Fig. 2). For a stable, symmetric (n^^ = n^), WIMP whose 
annihilation is s-wave dominated, mmin = 50 (10) GeV for / = I (0.2) if the particle is identical to its antiparticle 
(S: X = x); and rumin = 100 (20) GeV for / = 1 (0.2) if the particle differs from its antiparticle (NS: x ^ x)- Any 
further reduction in (CTti)cMB from future CMB experiments would increase mmin- For example, if the cosmic variance 
limit [13], a factor of four below the current Planck result [18] were reached, the lower bound on the DM mass would 
increase by a factor of four compared to the current constraints. 

These CMB constraints on WIMP DM candidates change dramatically if the WIMP is not a dark matter candidate, 
but only contributes a fraction of the DM mass density (fl^ < Udm, (o'i')x > (o’'1')dm)- As seen in §HIA and 


< 


(cru) 


CMB 


(cru) 


DM 


(9) 
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illustrated in Figs. 3-5, in this case there is no bound to the WIMP mass but, there is a lower bound to the 
annihilation cross section, ((t?;)!!!!!!, set by a combination of (tTu)DM and (crvjcMB- Although all masses are allowed, 
for > TOmin, (crw)min = (crw)DM, while for < minin, (av)min = (((^v)dm)^/(o-v)cmb > (crw)DM- The lower 
bound to the annihilation cross section corresponds to an upper bound to the relic mass density, < flmax, where, 
up to logarithmic corrections, ~ (o'w)min/(o'u)^ (see the right hand panel of Fig. 5). 

To illustrate the potential importance of the discussion here, consider the following application. There is interest 
in constraining the electric charge of a fractionally charged particle (q = Qe). A search for millicharged particles 
at STAG [25], led to upper bounds on Q for particles with masses in the range from 100 keV to 100 MeV. For 
WQ ^ 100 MeV < m^, annihilations can only lead to pairs or to photons. For mg ^ rrie and Q ^ 1, annihilation 
to pairs dominates, with an annihilation rate factor, {crv)Q « 2 x /itiq (cm^s“^), where mg is measured 

in GeV. At a mass of 100 MeV, the SLAG experiment set an upper bound to the charge of Q < 5.8 x 10“"^. Saturating 
this maximum charge, leads to an upper bound to the annihilation cross section, {av)g < 6.7 x 10“^® (cm^ s“^). 
Gomparing this upper bound to the GMB lower bound to {(Tv)g > (cru)i„in = 1-8 x 10“^^ (at rriy^ = 0.1 GeV) 
shown in Fig. 2, reveals a strong inconsistency. The GMB constraint on the annihilation cross section, in combination 
with the SLAG upper bound to the electric charge, eliminates the possibility of a 100 MeV millicharged particle. 
Indeed, extending the results for (avjmin presented here to lower masses and repeating this comparison, rules out any 
fractionally charged particle in the mass range from 100 keV to 100 MeV, complementing the constraints on higher 
mass fractionally charged particles (~ few GeV) presented in [26]. 

Before concluding, it is worth commenting on a recent preprint [27], in which Blum, Gui, and Kamionkowski also 
relax the assumption that the WIMP is a dark matter candidate, allowing for < Udm- The authors consider 
several observational consequences of late time annihilations, including the effect on the GMB. However, since they 
restrict their attention to larger WIMP masses than those investigated here, > 100 GeV > mmin, they do not 
consider the effects of the GMB constraint from late time annihilations on the lower bound (or not) of the WIMP 
mass and the lower bound to the annihilation cross section discussed here. As noted in [27] and here, the results 
for < Hdm have consequences for predictions of the expected gamma ray flux from late time annihilations in the 
present Universe in, e.ff., the Galaxy, dwarf galaxies, or the intergalactic medium in clusters of galaxies. In calculating 
the expected gamma ray fluxes, it is almost always assumed that the WIMP accounts for all of the dark matter in 
the astronomical target of interest (i.e., = tt-dm)- If, however, < udm, the predicted flux must be rescaled 

(reduced) by a factor of (Hj^/Hdm)^ < 1- In addition, the annihilation cross section must also be rescaled (increased) 
by a factor of (cru)x/(o’u)DM ~ Udm/^x ^ The overall effect is to reduce the expected gamma ray flux by a factor 
of (Oj^/Hdm) or, accounting for the logarithmic correction, by a factor of (Oj^/Hdm)]! — log (H;^/Hdm)/ 8]- The upper 
bound to this flux ratio, Umax/llDM, is shown in the right hand panel of Fig. 5. 

The analysis presented here has avoided the case of an asymmetry between particles and antiparticles (n^ 7 ^ n^) 
since for this case the relic abundance of the dominant particle (antiparticle) depends on the adopted asymmetry 
and not directly on the annihilation cross section. In the presence of an asymmetry, since the relic abundance of the 
subdominant antiparticle (particle) is suppressed by continued annihilations after the dominant particle has frozen 
out (n^ <C n^), the energy injection from the late time annihilations (pann cx n^n^) is also suppressed. As a result, 
it is not unlikely that an asymmetric WIMP could account for all of the dark matter, while being immune to the 
GMB constraint on late time annihilations (e.g., the lower bound to the WIMP mass could be much smaller than the 
value(s) of rumin derived here). For recent discussions of asymmetric dark matter, see, e.g., [28-30]. 

In summary, it has been shown here that if a WIMP %, a thermal relic, accounts for all of the dark matter 
(U^ = Udm), the GMB constraint on electromagnetic energy injection from late time WIMP annihilations sets a 
lower bound to its mass, > Wmin- For the current GMB data [18], this lower bound is m^m = 50 (100) GeV, for 
S (NS) particles and an electromagnetic energy efficiency factor / = 1 (for mmin ^ 10 GeV, mmin scales linearly with 
/ < 1). However, if the WIMP only accounts for a part of the dark matter, there is no limit to its mass, but there 
is a lower bound to its annihilation cross section, (avjann > (o'u)min, that increases inversely with the WIMP mass 
since, (cru)min = ((o-u)dm)^/(o'?;)cmb oc f/m^. 
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